Role of the transverse field in inverse freezing in the fermionic Ising spin-glass model 
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We investigate the inverse freezing in the fermionic Ising spin-glass (FISG) model in a transverse 
field T. The grand canonical potential is calculated in the static approximation, replica symmetry 
and one-step replica symmetry breaking Parisi scheme. It is argued that the average occupation per 
site n is strongly affected by V. As consequence, the boundary phase is modified and, therefore, the 
reentrance associated with the inverse freezing is modified too. 

PACS numbers: 



The inverse transitions (melting or freezing) , first pro- 
posed by Tammanni, are a class of a quite interesting 
phase transitions, utterly counterintuitive, in which the 
ordered phase has more entropy than the disordered one 2 . 
Despite this apparent unconventional thermodynamics, 
there is now a list plenty of physical systems in which 
this sort of transition appears (see Ref. [2| and refer- 
ences therein) including, interestingly, high temperature 
superconductors^ In that sense, the search for theoret- 
ical models which contain the necessary ingredients to 
produce such transitions has become a challenging issue 
as can be seen in Refs. Q and US However, much 
less considerations have been given to models in which 
quantum effects can also be taken into account. 

It should be noticed that there is an important differ- 
ence among some of the previous mentioned models. For 
example, in Refs. [2|, Schupper and Shnerb show that the 
Blume-Capel model can present inverse melting since the 
entropic advantage of the interacting state is introduced 
in the problem by imposing that the degeneracy param- 
eter r = k/l > 1, where it is assumed that ±1 spin states 
are /c-fold degenerated, while spin states are i-fold de- 
generated. In the Refs. and [j, the classical Ghatak- 
Sherringtonl (GS) model has also been proposed to be 
one of the simplest disordered models to present inverse 
freezing. Nevertheless, for the GS model, as remarked 
m Ref. H there is no need to enforce the entropic ad- 
vantage. Thus, the first order boundary of the spin-glass 
(SG) and/or paramagnetic (PM) transition naturally dis- 
plays a reentrance. That means it is possible to enter in 
the SG phase by heating from the PM one. That raises 
the following questions: is it possible to find other dis- 
ordered models with inverse freezing but no additional 
enforcement of the entropic advantage? Would it also be 
possible to incorporate quantum effects? 

It is now well known that the classical GS model has a 
very close relationship with the fermionic Ising spin-glass 
(FISG) 8 ' 9 model in the static approximation (SA) 10 . In 
the FISG mode l 11 ' 12 ! 13 , the spin operators are written by 
bilinear combination of fermionic creation and destruc- 
tion operators which act on a space with four eigenstates 
per site (|00), | f 0), |0 j), | U))- In particular, the ther- 



modynamics of both models can be exactly mapped by a 
relationship between the anisotropic constant D and the 
chemical potential /x of GS and FISG models^, respec- 
tively. Thus, we can expect strong resemblances between 
the phase diagrams of the two models; particularly, the 
presence of reentrance in the first order boundary must 
be emphasized (see, for instance, Ref. [Til ). Therefore, 
the FISG model can also be considered one of the mod- 
els which naturally presents inverse freezing. Actually, 
the FISG model has been intensively used to study the 
competition between the SG phase and, for example, su- 
perconductivity or Kondo effect (see Refs. [ID and [l6l 
and references therein). Moreover, the FISG has also 
been usecUi to study the effects of the transverse field T 
on the PM/SG boundary phase within SA. At the half- 
filling (fi = 0), it has been shown that the increase of T 
strongly changes the FISG thermodynamics since it tends 
to suppress the SG phase leading the freezing tempera- 
ture Tf to a quantum critical point (QCP) at r c . There- 
fore, the FISG model, besides allowing to treat charge 
and spin at the same level, can be a useful tool to study 
inverse freezing, particularly, when quantum effects can 
be included. 

The goal of the present paper is to investigate the in- 
verse freezing in the FISG model when spin flipping is 
induced by a transverse field T. It should be remarked 
that a set of FISG order parameters is composed not 
only by the usual nondiagonal SG order parameter q a p 
(a 7^ /?), but also by the diagonal one q aa , which is di- 
rectly related to the average occupation of fermions per 
site rai 3 -. Since T affects the behavior of the SG order 
parameters, we can assume that T could have strong in- 
fluence on n as well. Actually, the behavior of n is deter- 
mined by a more complicated dependence on T than the 
obvious one given by q aa . In that view, the spin flipping 
produced by T can also be considered as a mechanism 
to change the charge occupation and, hence, to modify 
the original phase boundary (when T — 0) in the T-fi 
plane. In this sense, we can probe the robustness of the 
inverse freezing when quantum effects are present, using 
what would be, in principle, a controlled external field. 
In this work, the partition function is obtained within 
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the Grassmann functional integral formalism^. The SG 
order parameters are calculated within the SA, in the 
replica symmetry (RS) and one-step replica symmetry 
breaking (1S-RSB) Parisi scheme. Previous calculations 
have showns that the position of the first order reen- 
trance has no significant difference when obtained in RS, 
1S-RSB or even in the full replica symmetry breaking 
(FRSB) Parisi scheme^. However, to be sure that the 
first order boundary is not affected by the RS instability 
when r is present, we also obtain the thermodynamics in 
the 1S-RSB. It should also be highlighted that the use of 
the SA in the present work can be justified since our main 
interest is to find the PM-SG phase boundaryi^ . The 
FISG model in the presence of transverse field is given 
by the Hamiltonian 



H 



(i) 



where the Jy coupling is a Gaussian random variable 
with mean zero and variance 16 J 2 . The spin operators 
in Eq. ^ are defined as Sf = \[n^ — and Sf — 

5[ C i T C H + C H C i T ] Witn "to = c\a C i<y ( CT =T) !)• We USe 

the procedure introduced in Ref. 1 - to obtain the grand 
canonical potential. Particularly, in the 1S-RSB Parisi 
scheme, the Grand Canonical Potential is written as 



/3fi = ^[(m - \)q\ - mql + q 2 } - fa 
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where 



K(z,v) = cosh{fa) + / D£cosh[y/A(z,v,g)], (3) 
withA(z,v,0 = [ph(z,v 1 0} 2 



((3T) 2 and 

h(z,v,Q = jV2(^/q^z + y/qx - q v + y/q-qiQ, 



(4) 



where Dx = dxe~ x ^ 2 /\Z2n (x = z, v or £). In Eqs. 
dU) and (g]), go an d qi are the 1S-RSB order parameters 
and q = q aa = (S^S^) is the diagonal replica spin-spin 
correlation. The parameters qo, q%, q, and m are given 
by the extreme condition of the grand canonical potential 
[Eq. [2] • The RS solution is recovered when qo = q\ (= q) 
and m = 0. In this case, the stability analysis of the 
RS solution is used in order to locate the tricritical point 
(Tt c ,Ht c ) (Ref. 21) as a function of T. Therefore, the 
condition for all eigenvalues of the Hessian matrix to be 
non-negative in the PM solution (q = 0) is 



q < 



for T/J > T tc /J 



- > MT,r,, W /,(T,r, g -)-4/(/w for T/J < TtJJ 
In Eq. ([5]), the tricritical temperature T tc is given by 



(5) 



T tc /J=^V2f Ttc (Ttc,r,q), 



(6) 



with ftf,(T, T,q) defined below, 



U(T,T,q) 





rjcj, cosh y/Aj, + sinh \fK$ 




J D£ (frl/Ai) cosh yTz; + (/3r)VAf ) 



where 



^={^ + 3(/3r) 2 [l-5(^ 2 /A )]}/A 2 
K0 =3(/3J) 2 [2h 2 - 1 + 5(/i|/A )] /A 5 / 2 , 



and 



Aa = h 2 



(7) 
(8) 

(9) 



In Eq. ©, f Ttc ( T tc T, q) is given from Eq. ©, w hen <ft = 
T tc and q = l/(y/2/3 J) which results in h Ttc = \A/2/3J£. 
Furthermore, there is no stable PM solution if 



fi< H at (T,r,q) for T/J>T tc /J, 



(10) 



M> H-(T,T,q) for T/J<T tc /J. 
In Eq. pop , the values of \i a t and /z_ are summarized by 
cosh -1 J D£[v<f, cosh y~K$ + sinh y/ A^] 



fi v (T,T,q) 
where 



V<t> = h 2 /(qA^) - 1, U0 = (/3r) 2 /(gAf ) 



(11) 



(12) 



with A^ and defined in Eq. ©. In Eq. (fTTj) . 

when ip ~ at, q = l/(v2/3J), and when </? = — , 
? = [/ (T,r,q) - V/0(T,r,g) -4/(/3J) 2 ]/2. 

In Eq. (jlOp . /i at defines the second order transition 
line T2/(/x) and /i_ gives the paramagnetic spinodal line 
[see Figs. l(a)-l(h)]. The fiT tc value is obtained by in- 
troducing T tc in equation for fj, at . 

The PM/SG phase boundary in the T/J-fi/J plane 
within RS, 1S-RSB Parisi scheme is plotted in Figs. 1(a)- 
1(d) for r = 0, 0.25 J, 0.5 J, and J, respectively. In Figs. 
l(e)-l(h) the position of (T tc , fj, tc ) (the tricritical point), 
the first order boundary phase and the spinodal lines are 
shown in detail. In particular, when T = 0, Fig. 1(a) 
displays the first order boundary confirming the existence 
of a reentrance in the FISG model. The turning on of 
the transverse field T [see Figs. l(b)-l(d)] produces a 
strong effect on the entire transition line. For instance, 
the second order part l2/(ju) is depressed. For /i = 
0, this behavior is reminiscent of the one already found 
in Ref. [17] in which the freezing temperature decreases 
toward a QCP when r — > r c . It should be remarked 
that the location of the tricritical point is also decreasing 
toward zero while T is enhanced. Nevertheless, the most 
important consequence can be seen in the behavior of 
the first order boundary phase Ti/(/i) [or equivalently 
fj,if(T)] and, in particular, in the reentrance [see Figs. 
l(e)-l(h)[. For T = 0.25J and 0.5J, T lf (fi) and spinodal 
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FIG. 1: Panels (a)-(d) show the phase diagrams /u/J versus T/J for several values of T / J, where T2/(/x) indicates the PM/SG 
second order phase transition, Tt c corresponds to the tricritical point, and Tif(fi) represents the behavior of the first order 
transition. The spinodal lines are also exhibited (full lines below Ttc)- Panels (e)-(h) show the first order boundary in detail. 
In these panels, Tif(fj) is presented for both the RS (dashed lines) and the 1S-RSB (pointed lines) solutions. They also exhibit 
the SG spinodal lines for the RS (dashed lines) and 1S-RSB (full lines) solutions. The dotted lines at very low temperatures 
are extrapolations to T = 0. 



lines are displaced in order to suppress the reentrance 
which is completely achieved when F = J. In Fig. 2, 
the grand canonical potential versus /i/J is plotted at 
T = 0.1J. This figure shows the displacement of first 
order boundary and spinodal points which illustrate the 



gradual suppression of the reentrance exhibited in Fig. 
1. In Figs. 1 and 2, results are shown within RS and 
1S-RSB schemes, which indicate that the location of the 
first order boundary is weakly dependent on the replica 
symmetry breaking scheme in agreement with Ref. U 
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FIG. 2: Grand canonical potential versus fi/J for several 
values of T/J for temperature T/J = 0.1. The dotted lines 
represent the PM canonical potential (Q pm ). The dashed and 
full lines represent the SG canonical potential (fl 3g ) in the RS 
and 1S-RSB solutions, respectively. The vertical lines are PM 
(left) and SG (right) spinodal lines. The full vertical line is the 
1S-RSB spinodal line, /ii/ indicates the first order boundary 
for T/J = 0.1. The labels SG pu and PM pu indicate the re- 
gions with only one spin-glass solution and one paramagnetic 
solution, respectively. 
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FIG. 3: Entropy versus T/J for several values of Y / J and \x/ J. 
The first order temperature is equal to Tif(fi) = 0.2 J for all 
figures. In panels (3a)-(3c), the left vertical lines are the SG 
spinodal lines and the right vertical lines T2/ are the PM/SG 
second order transitions. At low temperatures (T / J < 0.2), in 
(3a)-(3c), the dashed lines represent the RS spin-glass solution 
and the dotted lines indicate the 1S-RSB spin glass solution. 
In panel (3d), the vertical line is the PM spinodal line. The 
labels SGpu and PM pu indicate the regions with only one SG 
solution and one PM solution, respectively. 



There is some indication the T weaken even more such 
dependence. However, we can not be conclusive on this 
point due to numerical difficulties, in particular, at low 
temperature. 

In Fig. 3, the entropy versus temperature is show. The 
values of /i are adjusted to cross the point Tif(fi) — 0.2 J 
in the first order boundary transition. This procedure al- 
lows us to follow the entropy difference between SG and 
PM phases always at the same point of the first order 
boundary transition. For T = 0, we can see that the 
entropy of the PM phase is found below the SG one at 
the first order transition at ji — 0.92 J, which is expected 
from the existence of inverse freezing in the FISG model 
as previously discussed^ The increase in T produces the 
decrease in the entropy difference between SG and PM 
phases and the displacement of the first order boundary 
transition ^i/(0.2J) to larger values of fi. Simultane- 
ously, the pure PM phase (region PM pu in Fig. 3) is 
enlarged until the total disappearance of the pure SG 
phase (region SG pu in Fig. 3). Those effects are related 
to the suppression of the reentrance shown in Figs. 1 and 
2. 

The previous results clearly show that the spin flipping, 
due to the transverse field T, suppresses the reentrance 
in the first order boundary PM/SG transition and, there- 
fore, the inverse freezing. Actually, T depresses T2/(m) 
and the tricritical point which also implies that the first 



order boundary appears in a decreasing interval of tem- 
perature. It is known that, in the FISG model, the aver- 
age occupation of the nonmagnetic states exponentially 
decreases with the temperature^ Thus, the enhance of 
r would lead the phase transition to exist in a scenario 
where the nonmagnetic states are unimportant. That 
scenario is also consistent with the behavior of n as a 
function of T which has been studied by us for several 
isotherms. In that case, T tends to preserve the half- 
filling occupation even if \x increases. This effect becomes 
stronger when the temperature decreases. In that sense, 
the increase in T would redistribute charge in such way 
that the nonmagnetic states become gradually avoided 
at low temperature. There is no guarantee that the non- 
magnetic states in each site have been simply excluded 
by adjusting n. However, earlier results support such 
scenario. In Ref. [l?], the FISG model has been studied 
with an additional local restriction to get rid of the non- 
magnetic states. The obtained results have shown that, 
while r increases, the transition lines with the restriction 
and without restriction over the nonmagnetic sites at the 
half-filling become increasingly close. 

To conclude, in the present work, we have studied the 
role of spin flipping due to T in the inverse freezing of the 
FISG model. It should be remarked that FISG model 
presents inverse freezing when T = with no need of en- 
tropic advantage!^ As main result, it has been shown 
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that r destroys the reentrance in the PM/SG first or- 
der boundary and, thus, the inverse freezing. Our results 
suggest that T in the FISG model acts to redistribute the 
charge occupation, particularly, at low temperature. In 
that process, the nonmagnetic states become unimpor- 
tant for the phase transition. In that sense, T plays an 
opposite role concerning the inverse freezing as compared 



to the degeneracy parameter r of Ref. d. Although this 
results are restrict to the FISG model, we can speculate 
if the suppression of the inverse freezing by the increase 
in quantum effects can be a more general result. 
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